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Abstract. The extended nonlinear model has been appHed to construct neutron star 
matter equation of state. In the case of neutron star matter with non-zero strangeness 
the extension of the vector meson sector by the inclusion of nonlinear mixed terms 
results in the stiffening of the equation of state and accordingly in the higher value of 
the maximum neutron star mass. 



1. Introduction 

The analysis of the role of strangeness in nuclear structure in the aspect of multi-strange 
system is of great importance for both nuclear physics and for astrophysics. In the latter 
case understanding the properties of hyperon star is still a relevant item. At the core of 
a neutron star the matter density ranges from a few times the density of normal nuclear 
matter to about an order of a magnitude higher, at such densities exotic forms of matter 
such as hyperons are expected to emerge. The appearance of these additional degrees of 
freedom and their impact on a neutron star structure have been the subject of extensive 
studies HI- 

The description of a neutron star interior is modelled on the basis of the equation 
of state (EOS) of a dense nuclear system in a neutron rich environment. In general 
the description of nuclear matter is based on different models which can be grouped 
into phenomenological and microscopic. Additionally each one of them can be either 
relativistic or non-relativistic. In a microscopic approach the construction of the realistic 
model of nucleon-nucleon (NN) interaction can be inspired by the meson exchange theory 
of nuclear forces. The parameters within the model have to be adjusted to reproduce 
the experimental data for the deuteron properties and NN scattering phase shifts [5]. 
The solution of the nuclear many-body problem performed with the use of variational 
calculations for realistic NN interactions (for example for the Argonne vl4 or Urbana 
vl4 potentials) saturate at the density ~ 2 x po, where po denotes the saturation density 
[HI [7]. In order to obtain the correct description of nuclear matter properties, namely 
the saturation density, binding energy and compression modulus at the empirical values, 
a phenomenological three-nucleon interaction has to be introduced. Two-body forces, 
together with implemented three-body forces, help providing the correct saturation point 
of symmetric nuclear matter [8]. The nuclear matter EOS calculated with the use of 
the Brueckner-Hartree-Fock [9], [lO] approximation with the employed realistic two- 
nucleon interactions (the Bonn and Paris potentials) also does not correctly reproduce 
nuclear matter properties. Thus there are attempts to consider the nuclear interaction 
problem in a relativistic formalism. The relativistic version of the Brueckner-Hartree- 
Fock approximation - the Dirac-Brueckner-Hartree-Fock approach [TT] is also based on 
realistic NN interactions. The nuclear EOS obtained from the DBHF approach using 
the Bonn A potential is soft at moderate densities but become stiffer at higher densities 
[III Us]. At densities up to 2-3 times nuclear saturation density it is in agreement with 
constraints from heavy-ion collisions based on collective flow[T4l[l5] and kaon production 

m- 
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The relativistic approach to nuclear matter at the hadronic energy scale developed by 
Walecka is very successful in describing a variety of the ground-state properties of finite 
nuclei at, or near the valley of stability and in predicting properties of exotic nuclei with 
large neutron to proton excess. The standard Walecka model [17] comprises nucleons 
interacting through the exchange of simulating medium range attraction a mesons and 
u mesons responsible for the short range repulsion. Although this model properly 
described the saturation point and the data for finite nuclei it has been insufficient to 
properly describe the compression modulus of symmetric nuclear matter at saturation 
density and the proper density dependence in vector self-energy. The reproduced 
value of the incompressibility coefficient obtained in the original Walecka model gave 
large value of the order of ~ 500 MeV compared to the experimental results [18]. The 
nonlinear self-interactions of the scalar field (the cubic and quartic terms) have been 
added in order to get the acceptable value of the compression modulus [IHIEO]. Its 
estimation made on the basis of recent experimental data point to the range 230 ± 10 

MeV [211 [221 El EH . 

The quartic vector self-interaction softens the high density component of the EOS. This 
nonlinear term of the u vector meson has been used by Sugahara and Toki to construct 
TMl and TM2 models [25]. However, models which satisfactorily reproduce saturation 
properties of symmetric nuclear matter lead to considerable differences in the case when 
density and asymmetry dependence is included [26]. Isospin dependence of the strong 
interactions between nucleons infiuences both physical properties of nuclei and properties 
of infinite nuclear matter [27] . The latter case includes mainly the description of nuclear 
matter in high energy heavy-ion collisions and the properties of neutron star matter. 
Thus, the proper model of actual neutron star matter requires taking into consideration 
the effect of neutron-proton asymmetry. This in turn leads to the inclusion of the 
isovector meson p. The standard version of the p meson field introduction is of a minimal 
type without any nonlinearities. This case has been further enlarged by the nonlinear 
mixed isoscalar-isovector couplings which modify the density dependence of the p mean 
field and the energy symmetry. Such an extension of the neutron star model has been 
inspired by the paper [2H| in which the authors show the existence of a relationship 
between the neutron-rich skin of a heavy nucleus and the properties of a neutron star 
crust. 

The FSUGold model [21] which lead to considerably softer EOS. The meson sector of 
this model besides the linear terms of the scalar and vector fields includes nonlinear 
isoscalar meson self-interactions which soften the EOS of symmetric nuclear matter, 
and the mixed isoscalar-isovector coupling which alters the density dependence of the 
symmetry energy. Adding the mixed isoscalar-isovector meson interaction term the 
FSUGold model achieved acceptable results not only of the compression modulus for 
symmetric nuclear matter {K^ = 230 MeV) but the value of the neutron skin in ^°^Pb 
of Rn — Rp = 0.21 fm. The main astrophysical prediction of this model is connected 
with the value of the maximal neutron star mass which equals Mmax = 1.72M0. 
They are the nonlinear vector self-interactions that are discussed in this paper and the 



Nonlinear 



4 



hadronic SU{3) model which naturally includes nonlinear scalar and vector interaction 
terms is likely to be useful in the construction of models with nonzero strangeness and 
with more accurate description of asymmetric strangeness-rich neutron star matter [30] . 
The model considered has been extended to include a broad spectrum of nonlinear 
mixed vector meson couplings which stems from the very special form of the vector 
meson sector. The choice of these particular vector meson mixed interactions has been 
motivated by the chiral SU(3) model. The main effect of such an extension of the theory 
becomes evident when studying properties of neutron star matter, especially the form 
of the EOS. The equations of state for neutron star matter with hyperons considered in 
this paper have shown considerable stiffening for higher densities. Having obtained the 
equations of states the analysis of the maximum achievable neutron star mass for given 
class of models can be performed. Observational results limit the value of a neutron star 
mass and thereby put constraints on the EOS of high density nuclear matter. Recent 
observations of compact objects point to the existence of a high maximum neutron star 
mass [31]. The most spectacular result obtained with the Arecibo radio telescope for 
the neutron star-white dwarf binary system predicted the largest neutron star mass 
ever reported M{PSRJ0751 + 1807) = 2.1 ± 0.2MQ{la) [32]. However, this result 
has been corrected by Nice [131 IM]- The improved value of the orbital decay and the 
detection of the Shapiro delay lead to the new value I.26M0 but with the estimated 
errors 1.12 — 1.3OM0(1(t) and 0.98 — 1.53M0(2(t). Thus, this case can not be used as 
a constraint on the EOS. But there are another observations which indicate for high 
maximum mass of a neutron star. These are among others the value of the radius 
{R > 12 km)of the isolated neutron star RX J1856. 5-3754 [34] or neutron star mass in 
the low mass X-ray binary (LMXB) 4U 1636-536 estimated at the value of 2.0 ± O.IMq 
[35] . Another example of the LMXB is the neutron star in EXO 0748-676 which has 
been constrained by the detection of gravitational redshift of certain absorption lines. 
This combined with other observational data lead to individually estimated mass and 
radius of the star at the value of M > 2.10 ± O.28M0 and R > 13.8 ± 1.8 km [361137). 

2. Constituents of the model 

All calculations in this model have been done within the theoretical framework of quan- 
tum hadrodynamics (QHD) and the starting point is the nonlinear Walecka model which 
successfully describes the properties of nuclear matter and finite nuclei [17], [38]. Ow- 
ing to substantial asymmetry of neutron star matter models which include additional 
isospin carrying terms can be used for its remarkably complete description. Such models 
with Lagrangian functions supplemented by isospin dependent nonlinear, mixed vector 
meson couplings have been introduced and analyzed in papers by Piekarewicz et al. 
[2], [28], [41]. In the model considered there are couplings which relate the u and p 
vector mesons with the meson and thus link the asymmetry of the system with the 
strangeness content. The inclusion of a broad spectrum of mixed nonlinear vector me- 
son couplings provides the possibility of modifying the high density component of the 
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EOS. However, the presence of additional terms in the Lagrangian function requires the 
adjustment of new coupling constants. This has been done by fitting to the properties 
of nuclear matter. 

The presented model contains baryons and mesons as basic degrees of freedom. The 
baryon-meson interactions and meson-meson interactions are constructed on the basis 
of chiral SU(3) model [30]. 



2.1. Spin-O meson fields 

Additional hadronic states are produced in neutron star interiors in high density regime 
when the Fermi energy of nucleons exceeds the hyperon masses. The hadronic SU(3) 
|30| model which includes nonlinear scalar and vector interaction terms offers the 
possibility of constructing a strangeness rich neutron star model and of providing its 
detail description. The relevant degrees of freedom are hadrons - composite fields 
constructed from quarks. The problem considered is inspired by the chiral SU(3) model 
|30| in the nonlinear realization which infiuences transformational properties of quarks. 
The structure of hadrons is given in terms of their constituent quark fields q which can 
be split into left qi and right-handed qR parts {q = qi + qR)- They transform under 
5[/(3)l X SU{?,)r as 

(iL-*q'L = Lqi, qR^qR = RqR- (i) 

For the quarks of the nonlinear representation q the following relation can be written 

qi = uqi qR = v)qR (2) 

where qi and qR^x) are left- and right-handed components of the quark field q and 

u = u{j{x)) is given by 

«(vr(a;)) = gexp(^^vr"(x)A„^ (3) 

with pseudoscalar mesons tt^ considered as parameters of the symmetry transformation. 
They are identify with the octet of physical pseudoscalar fields [30] . 
In general the meson content of the model consists of spin-0 and spin-1 mesons. 
Introducing the matrix field $ enables the collective representation of the spin zero 
fields which in the matrix notation can be written as $ = S + ill. The S and 11 mesons 
can be transformed into nonlinearly transforming filds X and Y 

$ = S + m = u{X + iY)u, (4) 

where X is associated with the scalar nonet and Y with the pseudoscalar singlet which 
has to be added separately. The pseudoscalar mesons appear as the parameters of the 
symmetry transformation. The meson multiplets can be expanded in the basis of Gell- 
Mann matrices thus $ = ^T^^a where Ta = Xa are generators of U{3) and Xa (a=1...8) 
are the Gell-Mann matrices. It is convenient to introduce as a ninth matrix Aq = ^/ij'il 
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(/ is a unit matrix). The set {i\a}a=Q,...8 constitutes a basis of the Lie algebra u{3) of 
U{3). 

The scalar fields in the basis of U (3) generators are not mass eigenstates. If the ideal 
mixing between the octet and singlet states has been assumed then the octet state as 
and the singlet state cxo are related to the ideal mixing states a and C of a scalar nonet 
through the orthogonal transformation 



a \ I cos 6 sin \ ( 
(I I —sin 9 cos 9 J \ ctq 

Presenting the scalar multiplet as matrix the following form can be obtained 



(5) 



(-gO+g) 
«0 V2 



(6) 

In the process of spontaneous chiral symmetry breaking X acquires the vacuum 
expectation value (VEV) =< X >. As only components proportional to Aq and the 
hypercharge F ~ Ag are nonvanishing, < X > takes the form diag{< a >,< a >, < 
C >} and the following relations hold f-„ = \/ < a > and = (< ^ > + < C 
where and fx are the pseudoscalar decay constants. 

Making references to the Walecka model the following transformation should be done 
a- = (7+ < a > (7) 
C^a*+<C> 

where a and a* denotes fields in the Walecka model. 

In this paper the mean field approach serves as a method for solving the many body 
problem. In this approximation meson fields are separated into classical mean field 
values and quantum fluctuations, which are not included in the ground state. Thus, for 
the ground state of homogeneous infinite nuclear matter quantum fields operators are 
replaced by their classical expectation values Sq and Sq 

a ^ a + So (8) 
(7* = cr* + Sq- 



2.2. Spin-1 meson fields 

The spin-1 mesons are given by two octets of vector and axial vector fields. These fields 
similarly as in the case of spin zero mesons also can be written in a compact form 

W = 1{V,± A,) = ^ ± (9) 

' a=0 

where and corresponds to left and right-handed gauge fields, respectively and 
and Afj, denote the vector and axial vector nonets respectively. 

The physical isoscalar u and meson fields stem from the pure singlet and octet v^^ 
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isoscalar states and can be obtained assuming ideal mixing, namely that (j)^ is a pure ss 
state. This come down to the following relations 



0. = ^(a/2^J + <) 



(10) 



Accordingly the vector meson multiplet can be expressed explicitly in terms of physical 
fields as 

\ 



V2 



\ 



K*- 



*0 

K 



K*+ 

/ 



(11) 



2. 3. Baryon fields 

Nonlinearly transforming baryon fields can be written as 



(12) 



where and ^r are the left and right-handed parts of the baryon field in the linear 
representation. Baryon fields that enter the model are grouped into 3x3 traceless 
matrix B 

S+ 



B 



P 



n 



(13) 



3. The model 

The dynamics of the system has been described in terms of the Lagrangian function 
which in its most general form can be given as a sum of two basic parts representing 
meson and baryon sectors which are directly related by the term defining the baryon- 
meson interaction 

^ — ^kin'\' -\- -\- ^BM- (14) 

Additionally the kinetic term Ckin for both baryon and meson fields has been included. 



3.1. The meson sector 

The chiral SU(3) theory provides the basis for calculations made in this paper, however, 
the relativistic mean field approach to the description of the static, uniform nuclear 
matter leads to useful generalization about the model considered. In general, the meson 
sector includes contributions from spin zero and spin one mesons but in the mean field 
approximation vacuum expectation value of the pseudoscalar and axial meson fields 
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vanishes and these fields do not enter to the Lagrangian function. In the result meson 
field Lagrangian Cm embodies only the parts for scalar and vector mesons 

Cm = Cs + Cv (15) 
The scalar part of the Lagrangian function Cm includes the potential terms 

Cs = -h{l2f + k2h + 2k3l3, (16) 

which are given in the form of chiral invariants defined as 

Ji = Tr(X), h = Tr{X)\ h = detX. (17) 

The vector meson Lagrangian function presented in this model is the sum of a mass 
term and vector meson self-interaction terms up to fourth order in the fields 

Cv = \nlTr{V^V^) + Cyv (18) 



The second part of ( 1181) can be written in the form of invariants 

Cyy = ^ciTriV^V,))' + ^dTrliV^V.Y] + ^/[Tr(\/^)]^ (19) 

where is the vector meson multiplet. In order to split the mass degeneracy for the 
meson nonet the following chiral invariant has to be added 

CMv = \f^Tr[V,,V^''X% (20) 

This together with the kinetic energy term, which will be introduced in section 3.3, leads 
to the following terms for different vector mesons 



1. 1. 1 ,(T^ 

1. a\, ........ 1 



^[i - i^^-i^Wrf - ill - + ewi^^f (21) 



4[l-/^y](K.n'-4[l-/^fc'](^r)- 
As the coefficients are not equal unity the vector meson fields have to be renormalized 
by the factor Z^^{a) = 1 — /i(T^/2. The mass terms of the vector mesons differ from the 
mean mass niy by the renormalization factor and the following result can be obtained 

T^l = ^l = Z^{ao)ml; m\, = Zi^.(Co)"^y; = Z<^(Co)my (22) 

where the constants my /i and k are fixed to give the correct u and masses and 
o"o =< > |nfl=o cind Co =< C > |ns=07 dcuotcs the baryon number density. 
The baryonic part of the Lagrangian can be written as 

CB = Tr{BtYD^B), (23) 

where S is a 3 x 3 traceless hermitian matrix given by relation (fT3ll and denotes the 
covariant derivative of B which is defined as 

D^B = d^B + t[T^,B] (24) 

with defined as 

T/^ = ~2^^H9t^ + i9vlfi)u + u{df, + igvrf,)u^]. (25) 
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3.2. Baryon-meson interaction 

Using the notation introduced by [30] the very general SU(3) structure for the baryon- 
meson interaction can be written as 

Cbm = - V2gYiaw[BOBW]F + (1 - aw)[BOBW]D) (26) 

- gY^Tr{BOB)TrW, 

where W denotes general meson field and 

[BOBW]f := Tr{BOWB - BOBW) (27) 

[BOBW]d := Tr{BOWB + BOBW) - '^Tr{BOB)TrW. (28) 

Diff'erent forms of the presented above interaction terms result from differences in the 
Lorentz structure. 

For the nonlinear realization of chiral symmetry the antisymmetric (F-type) and 
symmetric (D-type) interaction terms of baryons not only with spin-1 but with spin-0 
mesons as well are allowed. 

In the case of baryon-scalar meson interaction W = X,0 = X, for vector meson 
W = V^,0 = 7^, for axial vector mesons W = A^,0 = and for pseudo-scalar 
mesons W = u^,0 = 7^,75. 

Baryon-spin zero meson interaction is indispensable for the construction of baryon mass 
terms. Masses of the whole baryon multiplet are generated spontaneously by the vacuum 
expectation values (VEV) of the non-strange and strange scalar fields. When the nucleon 
mass depends on the strange condensate < C > the parameters gf, g^ and as enable 
baryon masses to be fitted to their experimental values 

niMia, C) = "^0 - \gi{^as - 1)(V2C - ^) (29) 

"^a(o-, C) = 'Tio - '^giias - 1)(V2C - 0-) 

ms(a, C) = "^0 + \gl{^as - l)(v^C - ^) 

ms(a, = "^0 + \gl{'^as + 1)(V2C - f^), 

where mo is determined by two meson field condensates and is given by 

g V2< a> + <C> 
mo = g^ . (30) 

However, the assumption that = 1 and gf = y/Qgi leads to the model in which 
nucleon mass depends only on the nonstrange condensate < a >. In this case the 
coupling Constance between the baryons and the two scalar condensates are related 
to the additive quark model. Then nucleon mass depends only on the non-strange 
condensate < cr >, and only one coupling constant is needed to reproduce the correct 
value of the nucleon mass. To obtain the correct masses of the remaining baryons an 
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explicit symmetry breaking term has to be added. 

The explicit symmetry breaking term presented in the paper by [30] has been used 

jC-Am = -miTr{BB - BBS) - m^TriBSB) (31) 

where 5"^ = l/3(v^(A8)a — 6^) and the parameters g^, mi and m2 are used to determine 
baryon masses: 

rriN = - gNa < 0- > (32) 

1 _ 2 ^ - 

TTiH = - -gNa <cr > V2< ( > + "^i + 

6 ogNcr 

2 _ 1 /- - uii + 2m2 
rriA = - -gNa <(T > V2< C > H ^ 

6 ogNa <J 

2 _ 1 ^ - 
THE = - -gNa < a > V2< C > + mi. 

6 ogNa 

Considering the case when nucleon mass depends only on non-strange condensate and 
once again making references to the Walecka model the relation for the baryon masses 
can be expressed as 

TJiBia, C) = rriBicr, a*) = nxs - gsaCr - gBa*<y*, (33) 
where the terms gBaCr and gBa*cr* represent the modification of baryon masses due to 
the medium. 

The interaction of baryons with spin-1 mesons can be construct analogously with the 
baryon spin-0 meson interaction. For the case of pure F-type coupling (ay = 1) the 
assumption gf = V^gs (the strange vector field 0^ ~ S'j^s does not couple to nucleon) 
can be made. As it has been stated in the mean field approach the VEV of axial 
mesons are zero thus taking into account only vector mesons the Lagrangian describing 
baryon-vector meson interaction is given by 

Cbv = -V2gl {TrB^K, B] + TrBrB ■ TrV/:) (34) 
After insertion of the matrix ( ITTD to equation (i34ll and using all the facts concerning 
the construction of the baryon mass the following form of the Lagrangian function can 
be written 

Cms = Tr{B{tYD^ - mB{cT, 0)B) (35) 

where tub denotes the mass of the baryon octet in the chiral limit and the covariant 
derivative of B is given by 

D^B = d,B + t[r^,B], (36) 

with the connection F^. 

In the vector meson sector the couplings to the strange baryons determined from the 
symmetry relations read 

9ncu = (4ai/ - l)gs (37) 

9Aiv = 2 (Say - 2)^^ gA^ = — —{2av + 1)^^ 

gsiu = 2avg8 9i^<t> = -2V2(2ay - 1)51^ 
9elo = (2ay - l)gs 9r.^ = -2V2av5'^ 
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Taking ay = 1 the coupling constants related to the additive quark model can be 
obtained 

2 ^ V 9r.4, V2 , 

QKui — g^uj — gsui — -^9nuj — ^9s 9\<t> — 9t,4, — — -^9nuj- [oo) 

Similar symmetry relations can be obtained for the coupling constants in the scalar 
sector. These however are not used in the presented model. The couplings of baryons 
with the scalar mesons are determined from the experimentally estimated value of the 
A central potential. 



3.3. Kinetic terms 

The explicit form of the baryon and meson kinetic terms are given by 

Ck^n = tTr{B^^D^B) + ]^Tr{D^XD^^X) - ^Tr{V^,Vn (39) 
where denotes the covariant derivative 

D^ = d^ + i[T^,] (40) 

and is given by relation ( l25! l. The chirally invariant kinetic term for spin-1 mesons, 
in the case of vector mesons, is V^y = D^Vy — DyV^. The covariant derivative is 
given by relation fl40ll . 



4. The efFective model 



Vector mesons and thereby vector densities are the decisive factors that contribute to 
the EOS of dense matter in neutron star interiors. Thus, attention is focused on the 
construction of a model which includes wide spectrum of nonlinear couplings between 
vector meson fields. This allows one to perform a systematic analysis of their influence 
on the high density EOS. 

As has been stated in previous section theoretical description of strange hadronic matter 
requires the extension of the nonlinear Walecka model by the inclusion of baryons of the 
lowest SU(3) flavor octet. In order to describe the strongly attractive AA interaction two 
additional meson fields, the scalar meson /o(975) denoted as a* and the vector meson 
0(1020) have been introduced [39] . 

Thus, in the scalar meson sector besides non-strange a meson the hidden-strange a* 
scalar meson is included whereas, in the case of vector mesons cj, p and mesons are 
comprised. 

Summing-up the Lagrangian function for the system consists of a baryonic part which 
includes the full octet of baryons together with terms describing interaction of baryons 
with scalar and vector mesons and a mesonic part. The mesonic part contains also 
additional interactions between mesons which mathematically express themselves as 
supplementary, nonlinear terms in the Lagrangian function. Considering individual 
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constituents of the model which have been described in the previous paragraph the very 
general form of Lagrangian function can be written 

^BM = ^ ipBij^D^ipB - ^ ruBicr, a*)'ipj^ipB + ^m, (41) 



B B 

where baryon fields \E'^ = {ipN, ipA, ipT:, V's) ^ire composed of the following isomultiplets 



N 



is the covariant derivative of baryons which in terms oiu^.p'^^ and 0^ fields is given 

by 

Df, = + ig^B^^fi + ^g<t>B4>^l + igpBhB'r'^pl- (42) 
The meson part of the Lagrangian function 



Cm = \d^od^o + \d^a*d^a* - U,ff{a, a*, u^, p^, 0^) 
includes the field tensors Q^^u, ^fiu and i?^^ defined as 



(43) 

(44) 
(45) 



All meson interaction terms are collected in the potential function Uefficr, cr*,ujfj_, p'^^, (p^) 
which can be written as a sum of linear and nonlinear parts 

We//(a, (J*,UJ^, 0^) = Uiin{(T, a*,uj^, p^, 0^) + Unlio, d*, UJ^, p^, (j)^). (46) 

The linear scalar and vector meson part of the potential takes the form 



(47) 



imJ(pX'^)-im5(0,0^), 



whereas its nonlinear part is given by 

Uniia, cr*,u;^,p)^, 0^) 



Igsa' + \g,a' + h,c + 2d+ \m,r? + (48) 



+-^{c + d){plpn' + ^ciplp^M") + ^/(a;,0^)(0.0'^) + 

+\{c{pipn + ^d{pipn + 0(0,0^^) + ^/(0,0'^)(^.^^)) + 

+ i=/(0^a;'^)(a;,a;^) + \{c + d + D^u^u^^f. 
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The potential function We//(cr, o"*, cj^, p^, 0^) includes subsets of non-strange and 
strangeness rich matter and this can be written with the use of the following notation 

Uefficr, a\u^, p"^, 0^) = W,,5=o + ^,5^0, (49) 

where i denotes linear and non-linear parts of the potential (l46ll . If strangeness bearing 
components are not taking into account in the description of matter and only nucleons 
are considered, the Lagrangian density function constructed on the basis of the Walecka 
model can be retained which in the case of symmetric matter comprises contributions 
coming from nucleons, a and mesons. As neutron star matter is highly asymmetric 
one the inclusion of isovector-vector meson become indispensable. Finally making 
reference to the extended Walecka model with nonlinear scalar and vector meson self- 
interaction terms the Lagrangian function can be presented as follows 

C = ^^iil^D^, - (mAr(a) - gaN(T))i'N + ^df.ad^'a - 
+ I^PPIP"' - Ueff,s=o{(r, u;„ pi,) (50) 



2 



where 



N 




is the nucleon field, denotes the covariant derivative which now reduces to 

D^ = df, + ig^N^f, + igpNlmr^pl- (51) 

The nucleon mass is denoted by rriM^a) {N = n,p) whereas (z = a,uj,p) are masses 
assigned to the meson fields, i?^^ and fi^,^ again denotes the field tensors and are given 
by dm) and (l45ll. 

Collecting all the meson self-interactions and the mixed uo — p and a — p terms in the 
model considered a nonlinear part of the potential Ueff^s=o{<^,^fi:Plj) can be specified 

Unl,S=oi(T, U^, pD = - Jf^ + C3(^M^^)' (52) 



+^9'Ap>n' + iyv{g.NgpN?{u:,uj^Mpn + Ug.Ng,Nfa\plp^'^^ 



More detailed analysis of the vector meson influence on the high density EOS needs to 
consider different types of nonlinear vector meson couplings. 

The very general form of the vector meson Lagrangian (fT9l ) includes the potential which 
is defined by coupling constants that appeared in the vector meson Lagrangian (fTOll . 
Since the coupling constants of the mixed vector meson interactions are known very 
poorly, if ever, there are still important uncertainties in the analysis of their influence 
on the form of the EOS. In order to study the importance of individual parts of the 
potential it is necessary to find connections between the coupling constants c, d and / 
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and the couplings in the vector part of the nonlinear potential ?7„;(cr, cj^, p"). 



The results reveal the following relations 

c=^i3c,-3cs/3~2gXAv) (53) 

d=^i-cs + cs(3 + 2gXAv) (54) 

/ = C3/3. (55) 



In the paper by Serot et. al [40] the acceptable values of the parameters entering the 
very general nonlinear model have been introduced. In accordance with the estimations 
made in the paper |40| a new parameter (3 which takes the value in the range <0;1> can 
be introduced. This parameter makes possible to carry out more systematic analysis 
and to include wider class of models that have been constructed on the basis of the 
invariants presented in f fTOl ). The considered cases lead to models and results which 
already have been reported. Incorporating the relations (153II55I) the vector potential can 
be constructed 

= \cs {{p^pn'a + K^")') + (56) 
+ ioo.uj^) (^03(0^0^^) + ^-{g,g^yAy{2{p^y^) - (0^0^^)^ + 

+ (PlD (^C3(0^0^) - lf3<j>,r - lig.g.yAvM') + 
+ ^(0^0^)^ (^C3(l - ^f3) + 2{g^g,fAv^ 

+ Icsip'^^pniMm^ - /3) + ^f3cs{<Pur){<P,u;n + 

+ ^C3(0^0^)(^^^^) + -^PCs{cP,Cu'^){LU,CUn. 

Detailed analysis of the presented above potential can be made considering different 
cases. 

1. Nonstrange matter {S = 0). This case includes variety of models and leads directly 
to the well-known results presented in the literature. 

The potential (1561) reduces to the most general form for non-strange matter 

The analysis starts with the case P = and Ay = 0. 

Uv = \cs{p'',pn' + \cs{uj,ujn' (58) 

includes contributions from both isoscalar and isovector components. The last one 
is of special interest for the the asymmetric neutron star matter. 
The case with Ay 7^ 

l^v = Icsip^pn' + lcs{u;,u;n' + {g,g.f Ayip^pm^^^n- (59) 
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This form of the potential function incorporates additional vector meson mixed 
term. 

The case with /? = 1 and Ay 7^ 

Uv = \c,{uj,uj^f + {g,g^fKv{plpn{^,^n (60) 



This is the model studied by Piekarewicz et al. [28], [41] 
The case with (3 = 1 and Ay = 

Uv = icsKc^'^)' (61) 

The nonlinear Walecka model is retained in this situation. 

2. Matter with nonzero strangeness (S* 7^ 0). 
The case with /? = 

Uv = \cs{p;pn' + Icsip^pn^n + IcsiMn' m 
+ \cs{u^,un' - ^{g.g^rAvip^pnMn 

+ {gpg.fAviplpm^uUJn - \{gpg^fKy{ct^,r){^,^l 
The case with /? = and Ay = 

Uv = ^C3(pX'^)^ + ^C3(pX'^)(0.0^) 

+ ^c,i<p,^n' + lcs{<p,r)iu;.u;n (63) 



The case with /3 = 1 and Ay ^ 



+ ^cs{<P,u;n{oo,ujn + \c,{u;,u;n' - ]^{g,g^fKy{plpnM'') 



+ -^{gpg^fXv{ct>,rf 



+ {gpg^fKvP%uj,uj^) - ]^{gpg^fKv{(t),r){i^uUj'') (64) 
The case with (3 = 1 and Ay = 

Uv = ^cs{<p,rf + ^cs{<p,<pniru^u) + m 

3 11 

+ -C3(0^0'')(CJ,CJ") + -^C3(0^CJ^)(CU,CU'^) + -CsiuJ^^UJ^f. 

Summing up the vector meson potential in general can include a wide variety of 
nonlinear terms. The presence of nonzero strangeness led to very distinct division of the 
constructed models. This has been done with the use of the parameter /?. For nucleonic 
matter {S = 0) both cases (3 = and (3 = 1 have been analyzed. The latter one is of 
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special importance as it enables the comparison of the results with those obtained by 
Piekarewicz et al. |41| . However, in the case of nonzero strangeness important reduction 
has been done. Models constructed for the case (3 = 1 lead to very soft equation of state 
and consequently to very low value of the maximum neutron star masses. This is not 
interested according to results of recent observations which point towards larger masses. 
Such masses are especially concerned about the stiffening of the EOS at sufficiently high 
densities. This stiffening of the equation of state has become a big issue particularly for 
strangeness rich matter. Thus, further analysis of the influence of the nonlinear vector 
meson interactions for strangeness rich matter will concentrate on the case with f3 = 0. 



4.1. Mean Field Approximation 

The system considered, has been assumed to be isotropic, infinite matter in its ground 
state. To investigate properties of infinite nuclear matter, the mean field approximation 
has been adopted. The symmetries of infinite nuclear matter simplify the model to a 
great extent. The translational and rotational invariance claimed that the mean fields 
of all the vector fields vanish. Only the time-like components of the neutral vector 
mesons have a non-vanishing expectation value. Owing to parity conservation, the 
vacuum expectation value of pseudoscalar fields vanish (< tTq >= 0). Meson fields have 
been separated into classical mean field values and quantum fiuctuations, which are 
not included in the ground state. Thus, for the ground state of homogeneous infinite 
nuclear matter quantum fields operators are replaced by their classical expectation 
values. Hence, baryons move independently in the mean meson fields which generate 
themselves self-consistently by baryons. The resulting field equations for the mean field 
approximation have a reduced, simpler form 

ml-so + gssl + g^sl = ^g^Bmljj:j^S{meff,B,kF,B), (66) 

B 

^gujBUB, (67) 
B 

' f fp^O = Yl dpBhBUB , (68) 
B 

l'4 = Y9a'Bmljf^sS{meff,B), (69) 



B 

ml 

B 



B 



k = Yg^pBnB, (70) 



B 



where so,wo,ro,SQ and /o are the classical mean field values of the meson fields and 
rrii^eff are effective masses assigned to u, p and meson fields. The effective masses are 
given by the relations 

2 _ ™2 I Q„ „,,2 I OA ^ ^2^2 I n/"^ ^ ^ \ I „ „ ^2^ j-2 



+ Scsu^o' + 2ky{g^g,yrl + 2(-C3 - -Ky{g^gpY)f^ (71) 



'4 
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3 1 

mlf^p = ml + cgro + 2kv{gu^gpfwl + 2(-C3 - -Av{9u9pf)wofl (72) 



m,fr„ = m, 



+ Scsr^ + 2Kv{9ugpfwl + 4(-C3 - -Ay((?^(7p)2)w;^/o (73) 



The function S{meff^B, kp^B) is expressed with the use of an integral 



2Jb + 1 /-'^^ m.ff^B 



S{m,ff,B, kF,B) = / ,,, 7^'" fc^^fc, (74) 

where and I^b are the spin and isospin projection of baryon B, kF,B is the Fermi 
momentum of species hb denotes the baryon number density. The presence of the 
a* and meson fields provides new potential terms to the Dirac equation which now 
takes the form 

ih^df, - meff,B - gujBl^wo - gpBhBl^r^ro - g^Bl^h)i^B = (75) 

with nieff^B being the effective baryon mass generated by the baryon and scalar fields 
interaction and defined as 



meff,B = mB- {gaBSo + ga*BS*o). (76) 

In order to calculate the energy density and pressure of the nuclear matter the energy- 
momentum tensor T^^, which is given by the relation 

dC 

^^'^m^f'^'"''''^ ^^^^ 

have to be used. In equation (Fffll 0j denotes both boson and fermion fields. 
The energy density e equals < Too > whereas the pressure P is related to the statistical 
average of the trace of the spatial component Tij of the energy-momentum tensor. 
Calculations done for the considered model lead to the following explicit formula for 
the energy density and pressure: 

e = -ni^WQ + -C3W0 + -nipTQ + -m^/o + -m^.So + + (78) 
+^Av{gpgu)^wlrl + -^4 + U{so) + 3 Qcg + -^Av{gpgu.f] fo + 



+3 (jcs - ^Avigpg^^Y^ 



f2 2 



with es given by 



2 r'^F'B 

ee = ^ — / k^dk^/k"^ + {uib - g^BSo - g^^slY, 

_g ^ ^0 



(79) 



r> ^ 2,-'- 2|-'-/4|4\|-'-2f2 ^ 2 *2 /on\ 

P = 2^pro + 2^u.Wo + -csiwo + r^) + -m^f^ - -m^.s^ - (80) 
-U{so) + Avigpg^Ywy, + (^ca - ^Ayigpgu^Yj f^{wl + r^) + 



+ Qc3 + ^Ay(<7p^.)^)/o' + Pi3 
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The obtained form of the EOS determines the physical state and composition of matter 
at high densities. In order to construct the neutron star model through the entire density 
span it is necessary to add the EOS, characteristic for the inner and outer core, relevant 
to lower densities. Thus, a more complete and more realistic description of a neutron 
star requires taking into consideration not only the interior region of a neutron star, but 
also its remaining layers. In these calculations the composite EOS has been constructed 
by joining together the EOS of the neutron rich matter core region and neutron star 
crust. The inner crust is a region which spans from the neutron drip point to the inner 
boundary separating the solid crust from the homogeneous core [44] [45]. Since the 
density drops steeply near the surface of a neutron star, these layers do not contribute 
significantly to the total mass of a neutron star. The inner neutron rich region up to 
density p ~ 10^^ g cm~^ influences decisively the neutron star structure and evolution. 



5. The equilibrium conditions and composition of stellar matter. 

The ground state of a neutron star is thought to be the question of equilibrium 
dependence on the baryon and electric-charge conservation. Neutrons are the principal 
components of a neutron star when the density of matter is comparable to the nuclear 
density. For higher densities it is the equilibrium of the process 

p + e~ ^ n + z/g (82) 

which establishes the relation between chemical potentials 

Pp + fie = Pn + IJ-u,- (83) 

Thus, realistic neutron star models describe electrically neutral high density matter 
being in (3 equilibrium. The latter condition implies the presence of leptons. It is 
expressed by adding the Lagrangian of free leptons 

Y.iif{ird,-mf)i^f. (84) 

Neutrinos are neglected here since they leak out from a neutron star, whose energy 
diminishes at the same time. After electron chemical potential pe has reached the value 
equal to the muon mass, muons start to appear. Equilibrium with respect to the reaction 

e~ ^ pT + Ve + i^tj. (85) 

is assured when /i^ = pe (setting = /ip^ = 0). The appearance of muons reduces 
number of electrons and also affects the amount of the protons. 

Additional hadronic states are produced in neutron star interiors at sufficiently high 
densities when hyperon in-medium energy equals their chemical potential. The onset of 
hyperon formation depends on the attractive hyperon-nucleon interaction. The higher 
the density the more various hadronic species are expected to populate. They can 
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be formed both in leptonic and baryonic processes. In the latter the relevant strong 
interaction processes that establish the hadron population in neutron star matter e.g.: 

A + A^S + A^ (Q = 25MeV) (86) 

or 

S + A^^A + A^ (Q = SOMeV) (87) 

are Pauli blocked. Taking into account the energy released in these reactions (denoted 
as Q) it is likely that S hyperons do not appear in neutron star matter. The chemical 
potentials of neutron star components are related in such a way that the chemical 
equilibrium in stellar matter can be achieved. The requirement of charge neutrality and 
equilibrium under the week processes in the instance of strangeness rich matter 

B^^B2 + f + Uf B2 + f ^ B^ + Uf (88) 

leads to the following relations: 

E (>+ + ^/+) = E ("s- + ^/-) (89) 

i i 

where bi is the baryon number of particle i, qi is its charge, / stands for leptons f = e, fi, 
Bi denotes baryons and /Xj,^ = 0. The conditions mentioned above result in the relations 
between chemical potentials and constrain the species fractions in the stellar interior 
when taking into consideration the baryon octet and leptons included in this model: 

Aip = /iS+ = yUn - yUe /^A = /^SO = ^^=0 = IJ-n (90) 



6. Parameters 



Nuclear matter can be defined as an infinite system of nucleons with a fixed ratio of 
neutron to proton numbers and no Coulomb interaction. In general, the nuclear matter 
EOS, that is the energy per particle, of asymmetric infinite nuclear matter e(n6, fa) [46] 
defined as 

e{n,Ja) = - (91) 
rib 

is a function of two variables namely baryon number density Ub and the relative neutron 
excess fa (the asymmetry parameter) 

fa = (92) 

where rin and rip denote the neutron and proton number densities respectively. The sum 
nn + np = rib stands for the total baryon number density, £ in equation flOTj) denotes the 
total energy of the nuclear system. 

The properties of asymmetric nuclear matter can be studied with the use of the em- 
pirical parabolic approximation which allows one to expand the energy per particle of 
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Symmetric nuclear matter 




Symmetry energy 5*2 (nb) 


Co = e(no,0) 


j = S2{x = 0) 








= 27^3 {'-^] 


) 


K - r^i 



Table 1 . Coefficients that parametrize the behavior of the symmetric nuclear matter 
near saturation density. 



asymmetric nuclear matter in a Taylor series in fa 

eK, fa) = eK, 0) + 52K)/2 + S,{nb)f^ + ... (93) 

The factor fa makes the quartic 84(711,) term contribution negligible. 
Also the analysis performed with the use of realistic interactions indicates the dominant 
role of the 5*2 (n^) term not only in the vicinity of the saturation point but even at higher 
densities [9]. The expansion given in equation ([93l) enables the analysis of the function 
fa) in terms of the energy of symmetric nuclear matter e(nb, 0) and the symmetry 
energy S2{nb). Subsequently expanding e{nb, fa) around the equilibrium density Uq in 
a Taylor series in Ub, the following expressions for the two successive terms e{nb, 0) and 
5*2(^6) can be obtained: 



e{nb, 0) = e{no) + ]-K^x^ + \Qvx'^ + ... 
1 o 

S2{nb) = S2{no) + Lx + ]-KsymX^ + 



^X^' + ... 



(94) 
(95) 



where x denotes dimensionless parameter that characterizes the deviations of the density 
from its saturation value 

rib - no 



X 



3nn 



(96) 



Expressions (I94D and ((951) are parameterized by a set of coefficients: no, e(?T,o), 
J, Kgym, L and Qg which determine the behavior of the system near the saturation 
density. Particular coefficients are defined in Table[l]and evaluated at the point (no,0). 
This very point denotes the position of the state defined as the equilibrium state of 
symmetric nuclear matter e{nQ, 0) with minimum energy per nucleon and is characterized 
by the condition de{nb,0)/dnb = P(no,0) = 0. Thus, the linear term in the Taylor 
expansion (j94ll vanishes. 



Gathering altogether the terms of the expansions in n?, and in fa the approximated 
form of the EOS can be written as F46l 



eiribja) = e{no)+^{Ko+Ksymfa) 



Ub - no 
no 



+ 



Ub - no 



3 \ no 
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no 


av 




Qv 


TMl 


0.145 


-16.26 


281.16 


-258.295 


FSUGold 


0.148 


-16.27 


229.52 


-519.39 


RM 


0.148 


-16.0 


230.0 


-270.50 



Table 2. The parameters of the symmetric nuclear matter at saturation density. 





J 


L 


K 


Qsym 




TMl 


36.89 


110.79 


33.52 


-74.28 


-631.26 


FSUGold 


32.6 


60.44 


-51.51 


414.51 


-414.15 


RM 


32.0 


75.04 


-57.22 


343.98 


-507.46 



Table 3. The parameters of the symmetric nuclear matter at saturation density. 

Having obtained the EOS, each individual term that enters the formula (i97l ) can be 
calculated. 

According to the approximation presented by equation ( l93l l the symmetry energy can be 
calculated as the energy difference at a given density between symmetric {fa = 0) and 
pure neutron matter fa = 1. The density dependance of the symmetry energy around 
uq is determined by the parameters L and Kgym- Introducing the one-parameter fit to 
the low-density behavior of the symmetry energy 

where u = n/riQ and using this scaling properties the correlations between the density 
dependance of the symmetry energy and the neutron skin thickness can be estimated 
[471 [48l [49] . This dependance also allows one to determine the transition density pt 
between the crust and the core of a neutron star and to express it through the coefficients 
and Kgy^ in the following way 

The constraints on the value of 7t obtained from the intermediate-energy heavy-ion 
collisions provides a 7 value 7 ~ 0.69 — 1.05. Calculations performed in this paper are 
based on the standard TMl parameter set [25]. However, recent experimental results 
strongly indicate lower value of the symmetry energy coefficient and the compressibility 
coefficient of nuclear matter [48]. These lower values have been used to construct a 
parameter set (denoted by RM) which when compared with the TMl one, differs in 
the value of the scalar meson field mass. Also in the isovector sector the parameters 
gpN and Ay have been fitted to reproduce the symmetry energy coefficient at the value 
J = 32 MeV. The parameters and saturation properties of symmetric nuclear matter 
are collected in Table [2l 

The inclusion of the mixed nonlinear isoscalar-isovector coupling Ay provides the 
additional possibility of modifying the high density components of the symmetry energy 
and requires the adjustment of the QpN coupling constant to keep the same value 
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Ay 


ply 


J {MeV) 


L {MeV) 


K.„m i MeV) 


Kr (MeV) 


Ot ( fm ^) 


TMl foriff) 


0.0 


9.2644 


36.89 


110.79 


33.52 


-631.25 


0.1026 


TMl (nonl) 


0.0 


8.0642 


32.0 


96.12 


33.52 


-543.19 


0.1025 




0.008 


8.6567 


32.0 


85.13 


-37.98 


-548.77 


0.090 




0.01 


8.8264 


32.0 


82.38 


-50.93 


-545.23 


0.815 


RM 


0.0169 


9.234 


32.0 


75.04 


-57.22 


-507.46 


0.852 


FSUGold 


0.03 


11.767 


32.59 


60.44 


-51.51 


-414.15 


0.857 



Table 4. Saturation Properties of nuclear matter at saturation density obtained for 
nonlinear models 



of the symmetry energy at saturation. The remaining ground state properties are 
left unchanged. With these additional terms the expression for the symmetry energy 
coefficient Esymino) is now given by the equation 

^sym^^.) 8 mygl^ + 2Ky{g,^g^^)^wl 

H ^ (99) 

where where and niQ are the Fermi momentum and nucleon effective mass of 
symmetric nuclear matter at saturation. The first term in this equation coming from 
the explicit coupling between the nucleon isospin and the p meson whereas the second 
quantity is the relativistic kinetic energy contribution. The influence of the nonlinear 
couplings can also be considered in terms of effective u and p meson masses [50j which 
can be defined by the following relation 



m 



ml + 2Av{9pN9uNfrl, (100) 



m 



eff,P = ^l + '^^v{9pn9u^n)'^wI. (101) 

In this interpretation this is the p meson mass modification that influences the density 
dependence of the symmetry energy. The obtained form of the symmetry energy for the 
considered parameter sets are presented in Fig. [6l In general the nonlinearities soften 
the density dependance of the symmetry energy. For comparison the results of Akmal 
et al. have been included [51] 

Vector mesons-hyperon coupling constants are taken from the quark model and 
they are given by relations (l38l l. Whereas in the scalar sector the scalar coupling of the 
A and S hyperons requires constraining in order to reproduce the estimated values of the 
potential felt by a single A and a single E in the saturated nuclear matter. The analysis 
of the experimental data concerning the binding energies of A's bound in single particle 
orbitals in hypernuclei over an extensive range of mass numbers makes it possible to 
determine the potential depth of a single A in nuclear matter at the value of 

U[^^ - 27 - 30 MeV (102) 
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TMl (oi-ig) 
FSUGold 
RM 

TMl, A =0.0, J=32.0 
TMl, A^=0.008, J=32.0 
TMl, A^, =0.01, J=32.0 
APR 



50 - 




"0123456 

u=p/P„ 

Figure 1. The density dependence of the symmetry energy for nonHnear models. 
For comparison the results obtained for the ordinary TMl parameter set has been 
included. 

which corresponds to ~ 1/2 of the nucleon well depth f/]^''. There is still considerable 
uncertainty about the experimental status of E-nucleus potential. The calculations of S 
hypernuclei have been based on analysis of S~ atomic data. Phenomenological analysis 
of level shifts and widths in E~ atoms made by Batty et al. [52l [53] indicates that the 
S potential is attractive only at the nuclear surface, becoming repulsive for increasing 
density. The small attractive component of this potential is not sufficient to form bound 
S-hypernuclei. Also according to recent experimental data it has been established that 
the S nuclear interaction is strongly repulsive. 
The following values of the potentials have been used [55] 

U]^'^ = -28MeV, U^^^ + 30MeV, ui^^ = -ISMeV (103) 

for the determination of the g^A, QaT, and Qak coupling constants. In order to properly 
describe hyperon rich neutron star matter, the knowledge of the hyperon-hyperon 
interaction is indispensable. Data on AA hypernuclei are scarce. Observation of double- 
strange hypernuclei AA provide information about the A — A interaction. Several events 
have been identified which indicate an attractive AA interaction. The analysis of the 
data allows one to estimate the binding energies of ^^He, aaBg and aaB. The interaction 
between other type of hyperons are not known experimentally [541 [56] . The hyperon 
couplings to strange meson a* have been obtained from the following relations 

~ Uf> - 2U^^^ ~ 2U[^\ (104) 

In summery the coupling of hyperons to the strange meson a* has been limited by 
the estimated value of hyperon potential depths in hyperon matter this has also direct 
consequences for neutron star parameters. Recent experimental data [42] indicate a 
much weaker strength of hyperon-hyperon interaction. The currently obtained value of 
the uj^^ potential at the level of 5 MeV permits the existence of the additional parameter 
set which reproduces this weaker AA interaction. The parametrisation considered in 
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Y-Y interaction 


gaA 








weak 


6.17 


2.202 


5.41 


11.516 


strong 


6.17 


3.202 


7.018 


12.6 



Table 5. Strange scalar sector parameters 





Av = 


Av = 0.008 


Ay = 0.01 




9.2644 


10.207 


10.4828 



Table 6. The gpN parameters for different values of the parameter Ay 

this paper includes for comparison the strong and weak hyperon-hyperon couplings. 
The strong Y — Y interaction is related to the value of the potential f/|^^ = —20 MeV, 
whereas the weak one corresponds to uj^"^ = -5 MeV [42]. [43]. 

The inclusion of additional parameters in the isovector meson sector requires the 
adjustment of the QpN parameters. The new values of parameters are collected in Table 

m 

7. Results 

Having obtained the EOS that relates the energy density and pressure the corresponding 
solution of the Tolman-Oppenheimer-Volkoff (TOV) equations can be found and 
estimation of neutron star masses and radii become possible. The EOS and especially 
its high density limit has inevitable consequences for neutron star parameters. This 
manifests itself in a deep sensitivity of neutron star masses and radii to the stiffness of 
the EOS and allows one to study the influence of nonlinear vector meson interaction 
terms on neutron star properties. 

The integration of the TOV equations with a specific equation of state leads to 
the mass-radius relation and allows one to determine the value of the maximum mass 
which in a sense can give a measure of the impact of particular nonlinear couplings 
between vector mesons. In Fig. [2] the equations of state obtained for different cases 
of nonlinear potentials presented in this paper have been shown. Extreme, dashed 
curves represent results obtained for the standard TMl parameterisation, for the non- 
strange and strangeness rich matter respectively. The case when the matter includes only 
nucleons and leptons gives the stiffest EOS whereas the directly opposed result namely 
the softest EOS can be obtained for the standard TMl parameterisation extended by 
the inclusion of hyperons. The vector meson sector in this case comprises the quartic 
vector meson cu^ self-interaction term supplemented by the linear term for the hidden- 
strange meson 0^. The influence of the strength of hyperon-hyperon interaction is also 
illustrated by comparing equations of state calculated for the weak and strong Y — Y 
interaction. Other equations of state presented in this flgure aim to provide the analysis 
of the influence of additional vector meson nonlinear interaction. Relating this problem 
to the introduced scheme for nonlinear potentials, the consider cases Ay = 0.008 and 
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Figure 2. The EOS obtained for the nonlinear models. The limiting cases of the EOS 
indicate that the stiffest one has been obtained for the non-strange matter, whereas 
the softest represents the case of the standard TMl model extended by the inclusion 
of strange mesons which have been introduced in a minimal fashion. The remaining 
EOS for the nonlinear models with different values of the parameter Ay are located 
between these two curves. 

Ay = 0.01 are taking into account. The obtained results indicate for the strong tendency 
for stiffening of the EOS for the increasing value of the parameter Ay, approaching 
the limiting case for Ay = 0.01. Moreover, there exists evidence for diminishing the 
differences between the weak and strong Y — Y interaction for increasing value of the 
parameter Ay. Calculations performed for the value of the parameter Ay > 0.01 lead 
to acausal behavior of the equation of state at high density. 

The inclusion of nonlinear vector meson interactions has profound consequences for the 
structure of neutron stars and this can be deduced from Fig. [3] where the mass-radius 
relations for the obtained equations of state have been shown. Dotted curves depicts the 
results calculated for non-strange matter. Arrows marked by /5 = and P = I denotes 
the solutions obtained for ordinary TMl parameter set and for the TMl supplemented 
by the nonlinear uo — p coupling between the isoscalar and the isovector mesons. In the 
case of strangeness rich matter the results for different values of the parameter Ay has 
been included. The arrow marked Ay shows the influence of the nonlinear terms on 
neutron star parameters especially for radii. For comparison the mass-radius relation 
for the FSU Gold parameter set has been given. This figure depicts also constraints 
obtained from Vela and XTEJ 11739-285 data [58] . 

Fig [H depicts models obtained for the value of Ay = 0.008 with and without the 
nonhomogenous inner crust phase. This figure compares the solutions obtained for 
different structure of the outer layer of the neutron star. The presented results includes 
the homogenous neutron star model without the crust which is represented by dotted 
curve whereas dashed and solid lines depicts the mass-radius relations for the crust 
without and with the nonhomogenous inner crust, respectively. For both cases the 
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Figure 4. The mass-radius relation for the chosen models. It includes for comparison 
homogenous model without outer envelope. 



Ay 


Pmaxigcm-^) 






Ay = 


UA 


1.62 


12.92 


= 0.008 


27.4 


1.95 


10.53 


Ay = 0.01 


28.7 


2.161 


10.06 



Table 7. The value of the maximum mass configurations and the corresponding radii 
for different values of the parameter Ay 

critical density pt is given. A comparison of Fig. [3] and Fig. [4] leads to the conclusion 
that results obtained for higher value of the parameter Ay gives lower value of the 
critical density pt and as a consequence diminishes the nonhomogenous inner crust. 

The values of the maximum masses and the corresponding values of radii for the 
strangeness rich matter have been collected in Table [7l Results have been calculated for 
selected values of the parameter Ay. 
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Figure 5. The density profiles in the models representing the maximum mass 
configurations. 



The maximum mass increases with increasing value of the parameter Ay, given in 
the results neutron star models with masses exceeding 2Mq and with reduced radii. 
Thus, one can expect that solutions with nonlinear vector meson couplings lead to 
neutron star models with substantially greater density. The mass-radius diagrams have 
been obtained for variety of models presented in this paper. It includes the cases for 
non-strange and strange matter. Firstly solution for the pure TMl model only with 
quartic uj term is compared with the enlarged nonlinear TMl one which additionally 
comprises quartic p meson term. From these two mass-radius relations it is evident 
that the presence of the quartic p meson term influences mainly neutron star radii. 
Changing the value of the parameter Ay solutions with substantially reduced value of 
the transition density is obtained. This points to the conclusion that the crust-core 
boundary moves to lower density region leading to the models with reduced value of 
non-homogenous phase. This is confirmed in Fig. [5] which shows the density profile of 
neutron star matter for maximum mass configuration. Models with nonlinear vector 
meson couplings give as a result neutron stars with densities much more higher than 
that obtained for the case Ay = 0. 

The analysis of the density profiles indicates for the existence of different regions in 
neutron star interiors. These regions corresponds to the envelope and the core. In the 
core there is distinct part, with substantially increased density. This is connected with 
the appearance of strange matter. In the case of nonlinear models the inner core with 
nonzero strangeness spread through-out almost the whole star leading to very uniform 
neutron star model with considerably reduced outer parts. 

The appearance of hyperons follows from the chemical potential relations. It has been 
shown that the composition of hyperon star matter as well as the threshold density for 
hyperons, is altered when the strength of the hyperon-hyperon interaction is changed. 

However, the inclusion of nonlinear vector meson interactions also significantly 
modifies the chemical composition of the star. In order to get more complete 
understanding of the infiuence of nonlinear vector meson coupling on particular baryon 



Nonlinear 



28 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 

°0 5 10 15 20 

P'Po 

Figure 6. Relative concentrations of baryons for the model with Ay = 0, for the 
strong Y — Y interaction. 
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Figure 7. Relative concentrations of baryons for the model with nonlinear vector 
coupHng for Ay = 0.008, for the strong Y — Y interaction. 

and lepton concentration it is interesting to analyze this issue from different perspectives. 

Fig. [6] and [7] show relative fractions of particular baryon species YBi {i = 
n,p, A, S~, S~, S+ and S°) as a function of baryon number density Uf, for the strong 
value oi Y — Y interaction. Fig. [6] is constructed for the model without nonlinear 
vector meson coupling whereas Fig. [7] for the nonlinear model with Ay = 0.008. All 
calculations have been done under the assumption that the repulsive S interaction shifts 
the onset points of S hyperons to very high densities. 

A is the first strange baryon that emerges in hyperon star matter, it is followed by 
S~ and S° and S hyperons. These figures show differences in baryon fractions but 
these differences are more clearly visible in Fig. [8] and M presenting A and S hyperon 
concentrations for different models. 
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Figure 8. Concentrations of A hyperons for the cases discussed in the text, for different 
values of the parameter Ay . For comparison the results obtained for the weak hyperon- 
hyperon interaction is also included. The weak models are constructed for Ay = and 
Ay = 0.008. 
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Figure 9. Concentrations of S~ hyperons for the cases discussed in the text, for 
different values of the parameter Ay. For comparison the results obtained for the 
weak hyperon-hyperon interaction is also included. The weak models are constructed 
for Ay = and Ay = 0.008. 

Fig. [8] and M depict the analysis done for several models focusing on comparing the 
results for different values of the parameter Ay. Generally the appearance of nonlinear 
vector meson interaction terms lowers the concentration of hyperons in neutron star 
matter. Thus there is a substantial reduction of the A and H~ hyperon populations 
for the increasing value of Ay. In the case of S hyperon there is a density range for 
which the population of S hyperons vanishes. This results from the chemical equilibrium 
conditions which are set by relations among chemical potentials of the constituents of the 
neutron star matter. Chemical potentials depend on the effective baryon masses which 
have been substantially modified in the considered nonlinear models. For comparison the 
analysis of the infiuence of the hyperon-hyperon interaction strength has been included. 
Dashed lines in presented figures represent the concentrations of baryons for the weak 
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Figure 10. Baryon concentrations as a function of stellar radius for Ay — 0. 
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Figure 11. Baryon concentrations as a function of stellar radius for Ay — 0.008 
Y — Y interaction. 

In Fig. [10] and [TT] the analysis of the chemical composition of the star is depicted. 
These figures show profiles of particular baryon and lepton species as a function of 
the star radius. One can see that in the case of nonlinear models the hyperon core 
spreads through-out almost the whole interior of the star, however, it includes reduced 
population of hyperons. It is especially visible in the case os E~ hyperons. For the case 
Ay = 0.01, there is no E~ hyperons in the core. This fact ic strictly connected with the 
population of leptons, and especially influences the population of muons. 

From the presented flgures it is evident that lepton populations are altered not 
only by the change of strength of hyperon coupling constants but also by the isospin 
dependent nonlinearities which in turn determine the symmetry energy of the system. 
Thus a very special aspect of the existence of hyperons is the intrinsic deleptonization 
of neutron star matter. First the appearance of A hyperons stops the increase in the 
lepton population and additionally when negatively charged S~ hyperons emerge further 
deleptonization occurs. Thus the charge neutrality can be guaranteed with the reduced 
lepton contribution. On the contrary in the hyperon core of of the nonlinear models 
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Figure 12. Baryon concentrations as a function of stellar radius for Ay — 0.01 
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Figure 13. The asymmetry parameter fa and the strangeness fraction fs for models 
with Ay = as a function of the baryon number density. 

when negatively charged E~ hyperons disappear lepton concentrations is substantially 
enhanced, and large population of muons establishes. 

In Fig Us] and [H] the density dependence of the asymmetry parameter fa = 
{rin — Up) I {rih) which describes the relative neutron excess and the strangeness content 
fs, defined as (riA + 2n=- + 2n-^o)/nB are presented. Fig. [13] is constructed for the TM 
1 parameter set supplemented by the strange sector and for the case when the coupling 
Ay = 0. The second figure depicts the differences in the fa and fs behavior for two 
values of the Ay coupling. In the case that no nonlinearities are included the parameter 
fs takes the maximal value. Next the nonlinearities have been added with the result 
that the strangeness content of the system is reduced. The increase of the Ay parameter 
leads to even lower value of the strangeness fraction. The comparison of these two figure 
leads to the conclusion that the nonlinear model include matter more asymmetric but 
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Figure 14. The asymmetry parameter fa and the strangeness fraction for nonhnear 
models as a function of baryon number density. 
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Figure 15. The effective baryon masses for different value of the parameter Ay as a 
function of baryon number density. 
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Figure 16. The fraction of baryon and meson contribution in the total energy of the 
system. Calculations have been done for linear and nonlinear mdels. 



with lower value of the strangeness content. 
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Fig. [T5l presents changes of the effective baryon masses in dependance of the baryon 
density. It is common feature of the models that the effective baryon masses decreases 
with the increasing density. Very important aspects of the inclusion of nonlinear terms 
is connected with change of the effective baryon masses. In general the effective baryon 
masses are functions of sq and Sg. The numerical calculations show that the inclusion 
of nonlinear terms leads to the increase of the effective masses. Analysis has been done 
in comparison to the case when no nonlinearites are included. 

In Fig. [m the contributions of mesonic and baryonic parts to the total energy density 
is depicted. For the increasing density the nonlinear models lead to the situation when 
baryons and mesons give equal contributions to the energy of the whole system. 

8. Conclusions. 

Detailed knowledge about isospin asymmetric nuclear matter is of fundamental 
importance for understanding the structure of a neutron star, whose formation is 
preceded by the phenomenon of supernova explosion. In this paper a special class of the 
equations of state of asymmetric nuclear matter with non-zero strangeness have been 
analyzed in a systematic approach within the relativistic mean field model. The basic 
characteristic of the considered equations of state is the extended vector meson sector of 
the theory. This results in the appearance of various nonlinear vector meson couplings, 
among which there are terms which relate the strange and non-strange mesons. As 
a consequence strong connections between the asymmetry and strangeness fraction of 
the model have emerged. In order to construct neutron star models for the obtained 
equations of state the parameter sets which stem from the effective field theory and 
chiral SU(3) theory have been used. It has been shown that neutron star properties and 
through the properties also the neutron star structure are significantly affected not only 
by the presence of hyperons but also by the strength of hyperon-nucleon and hyperon- 
hyperon interactions. The results of the analysis performed for the nonlinear models 
have been compared with those obtained with the use of the standard TMl parameter 
set extended by nonlinear meson interaction terms, which have been added for detailed 
investigations of the high density symmetry energy. It has been shown that in the very 
nonlinear models the inclusion of hyperons does not soften the EOS; on the contrary 
it leads to its considerable stiffening. The consequences for neutron star parameters 
are straightforward and appear as a considerable growth of neutron star masses. Thus 
one of the inevitable conclusions is that in the case of nonlinear models the inclusion 
of hyperons does not results in the lowering of a neutron star mass. This is of special 
interest when considering pulsar data reported which indicate large neutron star masses 
and radii. This refers to measurements of the neutron star mass in pulsar-white dwarf 
system. The analysis performed clearly indicates that also the structure of a neutron 
star is changed in the case of nonlinear models. Stars become more uniform and more 
compact. The threshold for the appearance of hyperons is shifted to the very outer part 
of the neutron star core, but in general the hyperon fraction is reduced when compared 
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Figure 17. The dependance of the redshift z on the value of the parameter Ay 

with the linear models. The reduction of the hyperon population in neutron star matter 
is related to the lepton concentrations. The models show that for a sufficiently large 
value of the parameter Ay there are no S~ hyperons in the innermost part of the neutron 
star. Such particular models of the neutron star core reveal a hyperon reach layer with 
A and S~ hyperons and the central region in which S~ hyperons vanish. In the latter 
case the populations of leptons and especially muons get enhanced. 
The models analyzed include different types of nonlinear vector meson couplings. The 
additional coupling constants, which appear in these models are specified by the value 
of the parameter Ay which is constrained by the requirement of causality for neutron 
star matter. The value of Ay has profound infiuence on neutron star parameters. This 
is shown by Fig. [17] which depicts the dependance of the redshift z on the value of 
the parameter Ay. The presented models can be compare with those which do not 
include nonlinear vector meson interactions [59j, these solutions lead to very particular 
form of EOS This very particular form of the EOS generate different solutions of the 
Oppenheimer-Tolman-Volkoff equations. In the case of the cold neutron star model, 
apart from the ordinary neutron star branch, there exists an additional stable branch 
of solutions. 

The detection of redshifted O and Fe lines by XMM-Newton from the surface of 
the neutron star EXO 0748-676 [37] indicates for rather stiff EOS. 
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